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Abstract 

Using isospin relations, we predict the Standard-Model correlation between S^o Ks = (sin 2(3)^0 Ks and 
A^o Ks , the mixing- induced and direct CP asymmetries of B° — > ty°Ks- The calculation uses flavour 
SU(3) only to fix the isospin-3/2 amplitude through the B — > ty ± ty° branching ratio, and thus has 
a small irreducible theoretical error. It can reach percent level precision thanks to expected future 
lattice-QCD progress for the calculation of the relevant 5'L r (3)-breaking form-factor ratio, and serves 
as a benchmark for new-physics searches. We obtain an interesting picture in the A n o Ks ~S w o Ks plane, 
where the current experimental data show a discrepancy with the Standard Model, and comment on 
the direct CP asymmetries of B° — > ir~ K + and B + — > ty° K + . A modified electroweak penguin with 
a large new CP-violating phase can explain the discrepancy and allows us to accommodate also the 
corresponding data for other b — > s penguin-dominated decays. 
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Using isospin relations, we predict the Standard- Model correlation between S^o Ks = (sin 2[5) v o Ks 
and A n o Ks , the mixing-induced and direct CP asymmetries of B° — > ir°K$. The calculation uses 
flavour SU(3) only to fix the isospin-3/2 amplitude through the B —* tt 7r° branching ratio, 
and thus has a small irreducible theoretical error. It can reach percent level precision thanks to 
expected future lattice-QCD progress for the calculation of the relevant SU (3)-breaking form- factor 
ratio, and serves as a benchmark for new-physics searches. We obtain an interesting picture in the 
^-ttOk s ^S^ok s plane, where the current experimental data show a discrepancy with the Standard 
Model, and comment on the direct CP asymmetries of B° — > ir~K + and B + — » tv°K + . A modified 
electroweak penguin with a large new CP-violating phase can explain the discrepancy and allows us 
to accommodate also the corresponding data for other b — > s penguin-dominated decays. 
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Intriguing experimental results for observables of non- 
leptonic b — ► s decays [H have been receiving considerable 
attention for several years, where the 11 B — ► ttK puzzle" 
is an important example (see, e.g., @, 0, 0, 0, @, Q). 
The challenge is to disentangle possible signals of new 
physics (NP) from uncertainties that are related to strong 
interactions. In this context, a particularly interesting 
probe is offered by the time-dependent CP asymmetry in 



T(B°(t) -> ir°K s ) - T(B a (t) -> n°K s ) 
T(B°{t) -> ir°K S ) + T(B°{t) -> ir°K s ) 
= A„o Ks cos(AMd t) + S„o Ks sin(AAfd t) , 



(1) 



where S^o Ks arises from interference between mixing and 
decay, and A 7T ok s is the "direct" CP asymmetry. In the 
Standard Model (SM), we have - up to doubly Cabibbo- 
suppressed terms - the following expressions [8|: 



0, S^oxs = (sin 2/3) v o Kg * sin 2/3, (2) 



where (3 is one of the angles in the standard unitarity tri- 
angle (UT) of the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix. The current world average is [l[ 



(sin 2(3)^0 



A's 



0.58 ±0.17, 



(3) 



which should be compared with the "reference" value fol- 
lowing from B° — * J/ipKs and similar modes 



(sm2/3) J/ilKs = 0.681 ± 0.025. 



(4) 



The search for NP signals in the CP asymmetries of 
B° — > 7r /\s requires a reliable SM prediction of S n o Ks 
and/or A 7T o Ks . In this letter, we show that S^o^ can be 



calculated in the SM as a function of A 7r o Ks , with pro- 
jected irreducible theoretical errors at thepercent level. 
The starting point is the isospin relation [9(: 



V2A(B° -> tt°K ) + A(B° -» ir-K^ 
(f + 6)e*> + Awl -3A 3/2 ; 



(5) 



a similar relation holds for the CP-conjugate amplitudes, 
with A 3 /2 — * ^3/2 an( l 7 — > —7. Here T, C and P cw arc, 
respectively, the colour-allowed tree, colour-suppressed 
tree and electroweak penguin (EWP) contributions [lj|. 
The subscript of A 3 / 2 reminds us that the ttK final state 
has isospin / = 3/2, so that the individual QCD penguin 
contributions cancel in ([S])- S 7T o Ks can be written as 
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ir°K°) and An 



■sin(2/3-2^o Ks ), (6) 



If A 3 / 2 an( l A 3 /2 are known, 2<p w ox s 



B,rg(A 00 A* 00 ) 



can be fixed through ([5]), as shown in Fig. [TJ In order to 
determine j4 3 / 2 , we first rewrite the lower line of ([5]) as 



3^3/2 = -{T 



C)(< 



,»7 



(7) 



In the SM, the ratio qe luJ = -P ew /(f + C) is given by 



qe" 



-3 Co^) + CM 
2\ 2 R b Cx(ji) + Ci{n) 
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R q = 0.66 x — R q , (8) 
lib 



where A = \V US \ = 0.22, R b = 0.41±0.04cx \V ub /V cb \ is a 
UT side (value follows from flU), and the Cs are Wilson 
coefficients. If we assume exact SU (3) flavour symme- 
try and neglect penguin contractions, we have R q = 1 



2 




FIG. 1: The isospin relations §5§ in the complex plane. The 
magnitudes of the amplitudes, \Aij\ = \A(B — > K z tt :> )\ and 



IT, 



K*ir J )\, can be obtained from the corre- 



sponding branching ratios and direct CP asymmetries listed 
in Table HI while A 3 / 2 an d ^3/2 are fixed through ([8]) and ([9]). 



TABLE I: World averages of experimental data after 
ICHEP08 used in the numerical analyses (see also [l]]). 



Mode BR [1CT 6 ] 



A C T 



ScF 



ll|, |12[, while we shall use R q = 1 ± 0.3 for the numeri- 



cal analysis (results are robust with respect to the strong 
phase oj). Since qe lu factorizes at leading order (LO) in 
the 1 / mj expansion, R q can be well predicted using fac- 
torization techniques and future input from lattice QCD. 

SU (3) flavour symmetry allows us furthermore to fix 
\f + C\ through the b -> d decay B+ -> tt+tt [3]: 

\f + C\= R T+C \V US /V ud \ V2\A(B + 7r+7r°)|, (9) 

where the tiny EWP contributions to B + — > 7r + 7r° were 
neglected, but could be included using isospin [HI, [TBI] . 
We stress that ((9j) docs not rely on further dynam- 
ical assumptions. For the S'C/(3)-brcaking parameter 
Rt+c ~ Ik/ f-K we use the value 1.22 ± 0.2, where the 
error is quite conservative, as discussed below. 

Relations ([T])-© allow us to determine A 3 / 2 and A 3 / 2 , 
thereby fixing the two isospin triangles in Fig. [T] Since 
the triangles can be flipped around the A 3 / 2 and A 3 / 2 
sides, we encounter a fourfold ambiguity (not shown). 
Using ([6]), S 1T o Ks is determined as well. The correspond- 
ing prediction is shown in Fig. [2J where we keep A^o^ 
as a free parameter. For the implementation of this con- 
struction, we express the curves in Fig. [2] in parametric 
form |2j as functions of a strong phase 5 C , defined through 



r^e* 6 ' = {T + C)/P, 



(10) 



where P is the B° — > tt~ K + penguin amplitude [10j|. We 
find that no solutions exist for certain ranges of S c , sep- 
arating the full [0°,360°] range into two regions. They 
contain S c = 0° or 180° and correspond to the left and 
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FIG. 2: The SM constraints in the A 7r o Ks S n o Ks plane, as 
explained in the text. Left panel: contains <5 C ~ 0° (consistent 
with QCD), with S c — —60° (small circle), —30° (large circle), 
0° (star), 30° (large square), 60° (small square). Right panel: 
contains S c ~ 180° (not consistent with QCD), with S c = 
120° (small circle), 150° (large circle) , 180° (star), 210° (large 
square), 240° (small square). The shaded horizontal bands 
represent the value of (sin 2/3) j/^k s in Q. 
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FIG. 3: The constraints on r c e lS " that follow from the current 
data, as discussed in the text. Left panel: B — » ttK and B — > 
7T7r constraints (the symbols to label 8 C correspond to those in 
Fig- El- Right panel: B — > tth constraints for the BaBar and 
Belle data for A 7r + 7T - and the HFAG average. The solid and 
dotted lines refer to 1 a and 90% C.L. ranges, respectively. 



right panels of Fig. [21 respectively. As one circles the tra- 
jectory in either panel by changing 5 C , each value of this 
strong phase in the respective interval is attained twice. 
In order to illustrate this feature, we show - for central 
values of the input data/parameters - points correspond- 
ing to various choices of i5 c . The bands show the la 
variations obtained by adding in quadrature the errors 
due to all input data/parameters. Moreover, we assume 
7 = 65° ±10° [3G3- This angle will be determined with 
excellent accuracy thanks to CP violation measurements 
in pure tree B decays at the LHCb experiment (CERN). 

In order to resolve the fourfold ambiguity in Fig. [21 
we need further information on r c , S c : i) r c can be deter- 
mined if we fix \f+C\ through BR(B+ -> tt+tt ) (see ©) 
and |P| through BR{B+ -> tt+/v°) cx |F| 2 + . . . , where 
the dots represent negligible doubly Cabibbo-suppressed 
terms that are already strongly constrained by data [l8| . 
In the left panel of Fig. [3j the corresponding r c constraint 
is shown at the "charged" circle, ii) Using the SU(3) 
flavour symmetry and other plausible dynamical assump- 
tions 0] , a fit to all available B — > tttt data yields the tttt 



FIG. 4: The correlation in the A„o K „ -S„o K „ plane for a future FIG - 5: The SM correlation between A. 



■o K + 



-A, 



and 



benchmark scenario (narrow band) in comparison with the 
current situation (wider band), as explained in the text. 



curves. Since BaBar and Belle do not fully agree on the 
measurement of the direct CP asymmetry in B° — > tt + it~ 
[l| , we show in the right panel of Fig. [3] the correspond- 
ing allowed regions separately. We observe that the data 
imply S c ~ (0-30)°, in a gree ment with the heavy-quark 
expansion analyses in [3. Tl9j and [2(j, differing in their 
treatment of non-perturbative charm-penguin contribu- 
tions. Consequently, we can exclude the solutions shown 
in the right panel of Fig. [21 and are left with the twofold 
solution in the left panel. However, the lower band cor- 
responds to r c values of the "neutral" region in the left 
panel of Fig. [3] that are far off the right of the displayed 
region, drastically inconsistent both with the B — > tt-k 
data and with the heavy-quark limit. 

Consequently, we are left with the thin horizontal part 
of the upper band in the left panel of Fig. [2l which we 
show enlarged in Fig. [5] Using the experimental value 
for A 7T o Ks , we obtain the SM prediction 



S, 



°K S 



0.99 



+0.011 

-0.08 



hO.000 1 



hO.001 



cxp.-0.001|fl T+c -0.11|fl 5 



+0.001 
-0.07 



(11) 



which is about two standard deviations away from the 
experimental result in (|3|). It should be noted that (fTT|) 
depends on the input data collected in Table HI 

In Fig. [H we show the future theory error benchmark 
for the SM constraint in the A 7T o Ks ~S 1T o Ks plane. Both 
R q ([8]) and Rt+c ® factorize at LO in the 1 /m& expan- 
sion, and can be well predicted using input from lattice 
QCD. It should be stressed that "charming penguins" 
do not enter these ratios. As a working tool we use the 
approach of Ref . [1, [lj^ (BBNS) , but similar conclusions 
can be reached using Ref. [20j (where also derivatives 
of form factors would be needed). The key parame- 
ter is R q , which dominates the current theoretical error 
(fTTjl . Its uncertainty is governed by the S't/(3)-brcaking 
form-factor ratio & K = F B ^ K (0)/F B ^(0). If we as- 
sume = 1.2(1 ± 0.03), i.e. a 20% determination of 
the S'(7(3)-breaking corrections, as an optimistic - but 
achievable - goal for lattice QCD, we obtain the BBNS 



result R a 



(0.908±^)e l (°=;> o 



the present value R q 
larly, wc find Rt+c = 



(1.02 



to be compared with 

n+ 1 \ 



+0.27n l (0±i)° 



i 90+O.O2 
1 - zo -0.03' 



-0.227 



21|. Si 



A^o Ks for central values of inputs, with hadronic parameters 
fixed as for Fig. [2] (solid) , or following from the sum rule for 
rate differences [24f | (dashed). The dependence on S c is as in 
Fig. [5] and is constrained to SM values (upper curve in Fig. 
[2k.). 



precision is very mild as the form-factor dependence es- 
sentially cancels out. Setting, moreover, the uncertainties 
of the experimental inputs to zero, while keeping central 
values fixed, we obtain a prediction of S x ok s with errors 
at the percent level, as shown in Fig. [4] Consequently, 
the irreducible theory error of our proposed method for 
predicting S 7T o Ks in the SM is much smaller than in cal- 
culations using only the l/m& expansion, and makes it 
promising for a future e + e~ super-B factory (for a re- 
view, see, e.g., Ref. (22[)- 

Before turning to the interpretation of the current ex- 
perimental data in terms of NP, let us briefly comment on 
the difference of direct CP asymmetries A n o K + —A^- K +, 
which recently received quite some attention as a possible 
sign of NP [23j. Fig. O shows the SM correlation between 
this difference and the CP asymmetry A„ox s , keeping 
A 7T - K + fixed. It depends on CP-averaged B — > irK 
branching ratios and 7, and becomes equivalent to the 
sum rule for rate differences [24| when neglecting higher 
orders in sublcading amplitudes. We see that current 
data (cross) can be accommodated in the SM within the 
error on A^o Ks , although hadronic amplitudes then de- 
viate from the 1/mb pattern (see also Ref. [?})■ It would 
be desirable to reduce this uncertainty in the future. 

Let us now consider a NP scenario, which allows us to 
resolve the discrepancy between ^ and ([IT]) . Following 

0, we assume that NP manifests itself effectively in the 
data as a modified EWP with a CP-violating NP phase </>, 

1. e. q — > qe 1 ^ in Here q can differ from the SM value 
in l[8]). Since S c is rather small, the impact of this type 
of NP on A n o Ks and A n a K + is suppressed. In Fig. [SI we 
show constraints on qe 1 ^ from two \ 2 fits, using only the 
B — > irK data or both the B — > ttK and B — > ttit data. 
The latter have a strong impact on the allowed region 
of qe 1 ^ @, 0]> yielding two almost degenerate minima, 



5+0.2 



where the increase of 



q = 1.3±0.4, <j> = (63+g U )° and? = 0.8 
We also show the 90% C.L. regions (dashed curves) that 
correspond to a future scenario, assuming the benchmark 
value of Rq used in Fig. [5] and ten-times more data, with 
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FIG. 6: Constraints on qe 1 '*. Left panel: x 2 fit, using only the 
B — > 7rA' data. Right panel: x 2 fit, using both the B — > 7rA' 
and B — > 7T7t data. The inner and outer regions correspond 
to la and 90% C.L., respectively, while the stars denote the 
minima of the fits. The 90% C.L. regions with 10 times more 
data lie inside the dotted lines (see also the text). 




FIG. 7: Mixing-induced CP asymmetries for a set of penguin- 
dominated B° decays as functions of qsm(<fi), with qcos(4>) 
fixed to 0.6. The vertical bars depict the experimental 1 a 
ranges [l[ . The 1 a range (vertical band) and best-fit values 
(dashed line) for q sin <f> from Fig. [6] are also shown. 



relation in the A„ok s —S v ok s plane can be predicted reli- 
ably in the SM, with small irreducible theoretical errors, 
and have shown that the resolution of the present discrep- 
ancy with the data can be achieved through a modified 
EWP sector, with a large CP-violating NP phase. 
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